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I apply commonly used regularization schemes to a multiloop calculation to examine the properties 
of the schemes at higher orders. 1 find complete consistency between the conventional dimensional 
regularization scheme and dimensional reduction, but 1 find that the four-dimensional helicity scheme 
produces incorrect results at next-to-next-to-leading order and singular results at next-to-next-to- 
next-to-leading order It is not, therefore, a unitary regularization scheme. 

I. INTRODUCTION 

Dimensional regularization HI is an elegant and efficient means of handling the divergences that arise 
in perturbation theory beyond the tree level. Among its many favorable qualities it respects gauge and 
Lorentz invariance and allows one to handle both ultraviolet and infrared divergences in the same manner. 
The application of dimensional regularization to different kinds of problems has led to the development of 
a variety of regularization schemes, which share the dimensional regularization of momentum integrals, but 
differ in their handling of external (or observed) states and of spin degrees of freedom. 

The original formulation of dimensional regularization |lj, known as the 't Hooft-Veltman (HV) scheme, 
specifies that observed states are to be treated as four-dimensional, while internal states are to be treated as 
Dm = 4 — 2e dimensional. That is, both their momenta and spin degrees of freedom were to be continued 
from four to Dm dimensions. It turns out that one has the freedom to choose the value of the trace of 
the Dirac unit matrix to take its canonical value of four, so fermions continue to have two spin degrees of 
freedom, even though their momenta are continued to dimensions. Internal gauge bosons, however, 
have Dm — 2 spin degrees of freedom (internal massive gauge bosons have Dm — 1 degrees of freedom). 

A shght variation on the HV scheme has come to be called conventional dimensional regularization 
(CDR) fT\. In this variation, all particles and momenta are taken to be D^ dimensional. This often turns out 
to be computationally more convenient, since one set of rules governs all interactions. This is particularly 
so when computing higher order corrections to theories subject to infrared sensitivities, like QCD. In the 
HV scheme, if two external states have infrared sensitive overlaps, they must be treated as internal, or Dm 
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dimensional states. In the CDR scheme, all states are already treated as D,„ dimensional, so there is no 
possibility of failing to properly account for infrared overlaps. 

A third variation, called dimensional reduction (DRED) ||3l, was devised for application to supersym- 
metric theories. In supersymmetry, it is essential that the number of bosonic degrees of freedom is exactly 
equal to the number of fermionic degrees of freedom. This requirement is violated in the HV and CDR 
schemes. In the DRED scheme, the continuation to D„, dimensions is taken as a compactification from 
four dimensions. Thus, while space-time is taken to be four-dimensional and particles have the standard 
number of degrees of freedom, momenta span a D„, dimensional vector space and momentum integrals are 
regularized dimensionally. 

A fourth variation, the four-dimensional helicity (FDH) scheme H |5l, was developed primarily for 
use in constructing one-loop amplitudes from unitarity cuts. The most efficient building blocks for such 
calculations are tree-level helicity amplitudes, which necessarily have two spin degrees of freedom for both 
fermions and gauge bosons. The FDH scheme resembles the DRED scheme in that it regularizes momentum 
integrals dimensionally while maintaining the spin degrees of freedom of a four-dimensional theory (and 
therefore appears to be a valid supersymmetric regularization scheme ||5l), but there are crucial differences, 
which I will discuss in detail. 

The fact that the HV scheme respects the unitarity of the S'-matrix was proven at its introduction HI . The 
arguments which establish the validity of the HV scheme carry over to the CDR scheme and establish that 
it too is a valid regularization scheme. After some initial confusion over the proper renormalization proce- 
dure 10481 for the DRED scheme, it was established that it too is a proper, unitary regularization scheme 
and that it is indeed equivalent to the CDR scheme |9|. The FDH scheme has never been subjected to 
such stringent examination. It has been used successfully in a number of landmark next-to-leading order 
(NLO) calculations, but it has never been established whether it is a proper, unitary regularization scheme, 
or merely a set of shortcuts that allow expert users to obtain correct results. 

In this paper, I will perform a well-known multiloop calculation in the various regularization schemes. 
I will show that while the HV and CDR scheme calculations yield the correct result and the DRED scheme 
calculation, while far more complicated is completely equivalent, the FDH scheme calculation yields incor- 
rect results which inevitably violate unitarity at sufficiently high order. A detailed comparison of the various 
calculations identifies the source of the unitarity violations in the FDH scheme. 

The plan of this paper is as follows: in section two, I will describe the test calculation to be performed 
and present the result to be obtained. In sections three, four and five, I will describe in detail the calcu- 
lation to next-to-next-to-leading order (NNLO) as it is performed in the CDR, DRED and FDH schemes, 
respectively. In section six, I present partial results at N^LO which solidify the conclusion that the CDR and 



3 



DRED schemes are equivalent and correct, but that the FDH scheme violates unitarity. In section seven, I 
will discuss my results and draw my conclusions. 

n. THE TEST ENVIRONMENT 

To test the regularization schemes, I will calculate two quantities: the massless nonsinglet contributions 

to 

1. the hadronic decay width of a fictitious neutral vector boson V, of mass My; 

2. the single photon approximation to the total hadronic annihilation cross section for an electron - 
positron pair. 

I will perform these calculations by means of the optical theorem, taking the imaginary part of the forward 
scattering amplitudes. In both cases, this means taking the imaginary part of the vacuum polarization tensor 
sandwiched between external states. Since the optical theorem is a direct consequence of the unitarity of 
the 5-matrix, any unitary regularization scheme must give the same result, once one expands in terms of 
a standard coupling. To avoid comphcations involving prescriptions for handhng 75 and the Levi-Civita 
tensor, I will take V to have only vectorhke couplings. In this way, the vacuum polarization tensor for the V 
boson will be identical to that of the off shell photon, up to coupling constants and so the QCD expansion 
of the two results will differ only by constant numerical factors. 

Each regularization scheme will start from the same four-dimensional Lagrangian, 

jSf = - {d^d\l -^-^)-g^'^n)Al~gf'''id^'A'"')AlAl - ^f^'f'^'A^^A'^A'lA'^ 

. (1) 
+ j £ (5,7^ - igtU" - igv QfV) ¥} - c^Dc^ + gr''^ {d^ c") A'^ , 
f 

where A"^ is the QCD gauge field, is the massive vector boson, is the quark field of flavor /, c" and 
c" are the Faddeev-Popov ghost fields, g is the QCD couphng, gy is the V gauge couphng and Qf represents 
the charge of the quark flavor / under the V symmetry. I will not be computing nontrivial corrections in gy, 
so there is no need to specify the V-self interaction parts of the Lagrangian. 




FIG. 1: Sample diagrams of one-, two- and three-loop contributions to the vacuum polarization of V. 
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The result to N^LO is well known ifTOUUll . 
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To obtain the hadronic decay width at LO, NLO and NNLO, I need to compute the QCD corrections to 
the vacuum polarization of the V (photon) at 1, 2 and 3 loops, respectively. Sample diagrams are shown in 
Fig. ^. 



A. Methods 



In each scheme, I will need to compute the vacuum polarization of V and the necessary coupling renor- 
malization constants. As a cross-check on the reliability of my calculational framework, I reproduce known 
results on the QCD j5 -functions and mass anomalous dimensions to three-loop order, as well as the three- 
loop QCD contributions to the j3 -function of V (where needed). 

In all calculations, I generate the contributing diagrams using QGRAF |[T5]| . The symbolic algebra 
program FORM |[T6l is used to implement the Feynman rules and perform algebraic manipulations to reduce 
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the result to a set of Feynman integrals to be performed and their coefficients. The set of Feynman integrals 
are then reduced to master integrals using the program REDUZE lITTl . Using the method of Ref. flSl , the 
vertex corrections can be expressed in terms of the same propagator integrals used to compute the vacuum 
polarization and wave function renormalizations. The complete set of master integrals at one, two and three 
loops are shown in Fig. Q. Most of the master integrals are trivial iterated-bubble diagrams and the others 
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FIG. 2: Master integrals for the evaluation of vacuum polarization at a) one loop, b) two loops and c) three loops. 

were evaluated long ago |[T9ll20l . As an additional cross-check, the integral reduction and evaluation is also 
performed using the program MINCER lETl l22ll . 

B. Notation 

The various schemes that I will consider span a variety of vector spaces, each with their own metric 
tensor. To establish some level of consistency, I will denote the metric tensor of classical four-dimensional 
space-time as T]^^; the metric tensor of the D„, dimensional vector space in which momentum integrals 
are regularized will be denoted as g^^; and the metric tensor of the largest vector space will be denoted 
gf^^. Where it does not vanish, the complement of g^^^ will be denoted as 5^^ = g^^ — g^^- Similarly, the 
Dirac matrices 7^ , will be denoted when they are strictly four-dimensional, 7^ when they span the 
dimensional space and in the space spanned by 5'^^. 

I will now present the details of the calculation in the CDR, DRED and FDH schemes. 
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in. CONVENTIONAL DIMENSIONAL REGULARIZATION 

In the CDR scheme, the calculation is quite straightforward. The Lagrangian and Feynman rules are just 
the same as for a four-dimensional calculation, except that the Dirac matrices and the metric tensor g^^^ 
have been extended to span a Dm dimensional vector space. That is. 

The Dirac trace, Tr [1] = 4, retains its standard normaUzation. 

Although Dm is given the representation = 4 — 2 e, the sign of e is not determined. If it is taken to 
be positive, so that D,„ < 4, then the Feynman integrals that one encounters are convergent under the rules 
of ultraviolet power counting. On the other hand, infrared power counting would prefer e < =^ D,„ > 4. 
In practice, the sign of £ does not matter and it can be used to regularize both infrared and ultraviolet 
divergences. Regardless of the sign of e, it is important that the vector space in which momenta take values 
is larger than the standard 3 + 1 dimensional space-time. This means that the standard four-dimensional 
metric tensor T]'^^ spans a smaller space than the Dm dimensional metric tensor, and the four-dimensional 
Dirac matrices y^'i'^'^ form a subset of the full y'^, 

e'8'^=8''. rnl = n''. r^^^ rjP = r^^P . (5) 

These considerations are of particular importance when considering chiral objects involving 75 and the 
Levi-Civita tensor, but will play a role in our discussion below. 

Because the Dirac trace is unchanged, fermions still have exactly two degrees of freedom in the CDR 

scheme. Gauge bosons, however, acquire extra spin degrees of freedom in the D,„ dimensional vector space. 
The spin sum over polarization vectors in a physical (axial) gauge takes the form 

-g^vIe*^A:,A)e^(fc,A)=g^v(g^^- \=Dm-2 = 2-2e, (6) 

where n is the axial gauge reference vector. For massive vector bosons, the spin sum becomes 



(7) 



A. Renormalization 



The renormalization constants in the CDR scheme are defined as 



^ AAA — ■^i'- AAA ^ V/ —^l^f-' — ^3 
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where Fate represents the vertex function involving fields a, b and c. 

Although we treat the quark fields as massless, we can compute the mass anomalous dimension by 
introducing a fictitious scalar particle <p and computing the jS-function of its Yukawa coupling to the quarks. 
The equivalence is clear from the standard model, where the Higgs Yukawa coupling and the fermion mass 
are proportional at leading electroweak order and must behave the same under QCD renormalization. For 
this purpose, I introduce one more renormalization constant, F^^^ = Zi^Fg^ip. One can introduce a wave 
function renormalization for 0, Z30, but it will not contribute because = 1 + (^(a^). Note also that I do 
not need to compute the QCD corrections to the j8 -function for ay, which will start at order ay because of 
the Ward Identity. 

In the MS scheme, the couplings renormalize as 

' ^<^-z|-zp3-ifz3 



The structure of the renormalization constants Z m and Z hs is determined entirely by their lowest order 
(1/e) poles, which in turn define the j3-functions. 
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The mass anomalous dimension, 

^ ' > m^sd^^ y 71 J 

is defined in terms of m, rather than m^, with the result that = ^j8^^. The results for jS^^ and 
through three loops are given in Appendix [A] 
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B. Vacuum polarization in the CDR scheme 

The imaginary part of the unrenormalized vacuum polarization tensor in the CDR scheme is 

,2 . 
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Upon renormalizing the QCD coupling according to Eq. (9 1, setting tty — )• ay ^ ^^4^ j > dropping terms 
of order {s),l obtain 
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In this way of performing the calculation, all of the QCD states that appear are internal states, so the HV 
scheme gives exactly the same result. 

C. Total Decay rate and annihilation cross section in the CDR scheme 



The decay rate and the annihilation cross section are determined by computing the imaginary part of the 
forward scattering amplitude. For the decay rate, this means attaching the polarization vector e'^{Q,X) and 
its conjugate e^(Q, A)* {Q^ = My) and averaging over the spins. 
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Notice that because the imaginary part of the vacuum polarization tensor is finite, it does not matter 
whether the spin sum is taken in D,„ = 4 — 2 e dimensions as in the CDR scheme or in four dimensions as 
in the HV scheme as the difference is of order e. The result is 
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in agreement with Eqs. (2]|3 



For the annihilation cross section CJ^+g-^ hadrons » one attaches fermion bihnears to each end of the vac- 
uum polarization tensor and averages over the spins. 
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Because this is a forward scattering amplitude, the spinor bilinears can be combined into a trace, 

\Y^{v{p,.^)\y^\u{p,-^'))(u{p,-^')\f\v{p,.^)) = \T^^^ 

XX' 

(18) 

where the last identification results from the fact that = + p^+, pg- ■ Q = pg+ ■ Q = The 
result is 
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again in agreement with Eqs. (|2]j3]). 

Thus, I have established that I can reproduce the known results in the CDR scheme through three-loop 
order, which is a strong check on my computational framework. 



IV. DIMENSIONAL REDUCTION 



In dimensional reduction, one starts from standard four-dimensional space-time and compactifies to a 
smaller vector space of dimension D,„ = 4 — 2£<4in which momenta take values. The particles in the 
spectrum, however, retain the spin degrees of freedom of four dimensions. That is, both fermions and gauge 
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bosons still have two degrees of freedom. This is by design, of course, since it is required by supersymmetry. 
All Dirac algebra can be treated as four-dimensional. However, now the four-dimensional metric tensor t]^^ 
spans a larger space than the D„ dimensional metric g'^^ that might arise from tensor momentum integrals, 

rn'^=g'p. (20) 

There is also a very serious consequence of the fact that the dimensional vector space is smaller 
than four-dimensional space-time. The Ward Identity only appUes to the dimensional vector space. 
This means that the 2 e spin degrees of freedom that are not protected by the Ward Identity must renor- 
malize differently than the 2 — 2 e degrees of freedom that are protected. In supersymmetric theories, the 
supersymmetry provides the missing Ward Identity which demands that the 2e spin degrees of freedom be 
treated as gauge bosons. In nonsupersymmetric theories, however, they must be considered to be distinct 
particles, with distinct couplings and renormalization properties. It is common to refer to these extra degrees 
of freedom as "e-scalars" or as "evanescent" degrees of freedom. 

Once the evanescent degrees of freedom (which I will label A"^, to distinguish them from the gluons, 
A"/^) are recognized as independent particles, it is apparent that their couplings are also independent, not 
only of the QCD coupling, but of one another. That is, the coupUng ge of the evanescent gluons to the 
quarks is not only distinct from g, the coupling of QCD, but is also distinct from A;, the quartic couphngs 
of the evanescent gluons to themselves. (The quartic gauge coupUng of QCD sphts into three independent 
quartic couplings of the evanescent gluons.) Note that the massive vector boson also has evanescent 
degrees of freedom, Vf, which couple to quarks with strength gve- 

Thus, the Lagrangian in the DRED scheme becomes: 

+ i^U^ - iStf/" - igvQfV) ¥j - c^a + gf' {d^, c«) A"^ 

f 

+ Ia^^D Af - gf''{d^Af)A^^AU + ^^f^ffA'^AfAiAi^ - \ Y^^iH^''^ Af Af A%Ai, 

f 

(21) 

As mentioned above, the quartic coupling of the evanescent gluons splits into three terms, which mix under 
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renormalization. One can choose the tensors Hf'^'^'^ to be ||23]| 
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Although the quartic coupUngs enter the jS -functions and anomalous dimension at three loops and are es- 
sential to the renormalization program, they do not explicitly contribute to the calculation at hand. 

Now that the correct spectrum has been identified, one must carefully consider the renormalization 
program. The naive application of the principle of minimal subtraction leads to the violation of unitarity IH. 
Because the contributions of evanescent states and couplings to scattering amplitudes are weighted by a 
factor e, the leading one-loop contribution is finite and therefore not subtracted. As one proceeds to higher 
orders, there is a mismatch among the counterterms such that the renormalization program fails to remove 
all of the ultraviolet singularities. 

A successful renormalization program for the DRED scheme lIH |9l applies the principle of minimal 
subtraction to the evanescent Green functions (that is. Green functions with external evanescent states) 
themselves. At each order, the renormalization scheme renders the evanescent Green functions finite. Since 
evanescent Green functions enter into the scattering amplitudes of physical particles at order e and they are 
rendered finite by renormalization, they never contribute to physical scattering amplitudes. 

The evanescent coupling still contributes to Green functions with only physical external states, but the 
contribution is rendered finite by the prescribed renormalization program ||8l|9l|23l|24l. Because the evanes- 
cent coupling, ttg renormalizes differently than the gauge coupling a,, the two cannot be identified, even at 
the end of the calculation. One can choose a renormalization point where the two coincide, but they evolve 
differently under renormalization group transformations and their values will diverge as one moves away 
from the renormalization point. 

Still, the evanescent coupling is essentially a fictitious quantity and one finds that if one computes a 
physical quantity in the DRED scheme and then converts the running couplings of the DRED scheme to 
those of a scheme such as CDR that has no evanescent couplings, the factors of drop out ||23ll24ll . 
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A. Renormalization 



The renormalization constants in the DRED scheme are defined as 
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In addition, I again introduce the fictitious scalar that allows me to compute the mass anomalous dimension 
for massless quarks. Note that while the Ward Identity protects ay from leading QCD corrections, it does 
not protect Ove- That is why I need to introduce renormahzation constants for the vertex and wave-function 
and why I need to compute the j8 -function of Uve- 

In the DR scheme (modified minimal subtraction in the DRED scheme), the couplings renormalize as 
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and the jS -functions are given by 
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Through three-loop order, the T], do not contribute to the QCD jS -function, nor to the vacuum 
polarization of V (or Ve). To three-loop order, I find agreement with known results Il23ll24ll and derive new 
results for the jS -function of ay^. The coefficients of the j3 -functions and anomalous dimensions are given 
in Appendix [B| 



By comparing and y in Eqs. (B4 



B5l, we see that the term "£-scalar" is a misnomer. If the 
evanescent part of V were a true scalar, its j3 -function would coincide (but for a factor of 2) with the mass 
anomalous dimension. The pure af^ terms do coincide, because there is no nonvanishing contraction of 
the Lorentz indices of the evanescent V and those of the gluons. Because there are contractions between the 
Lorentz indices of the evanescent V and those of the evanescent gluons, however, terms involving do 
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not agree. 

Calculations in the DRED scheme naturally produce results in terms of while the standard result 
has been expressed in terms of af^^. One can always convert one renormaUzed coupling to another. The 
rule for converting af^ — )- a^^, derived in Refs. ||24l|25l, is 
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When the result is expressed in terms of q^^, all af^ terms drop out. 



B. Vacuum polarization in the DRED scheme 



In the DRED scheme, there are two independent transverse vacuum polarization tensors, 
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Upon renormalization according to Eq. (24 1 and expanding in terms of according to Eq. (26 1, 1 find 



16 



that 



DR 



/ 



+ 



a 



DR 



n I 4 
121 11 



a 



DR 



n 



I rdr. M 

|P20 1117^ 



,3 / 121 11 \ /111 



1 + 



a; 



MS 



/3o^^ In 



+ 



MS 



1Z 



123 lU , ^ 



11 1 „ 
16 + 2^^ 



3[nB(e)| 



DREDi 



a, 



MS 



C. Total Decay rate and annihilation cross section in the DRED scheme 
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As in the CDR scheme, the decay rate and annihilation cross section are determined from the imaginary 
part of the forward scattering amplitude. 
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For the annihilation cross section Og+e-^ hadrons^ one attaches fermion biUnears to each end of the vac- 
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uum polarization tensor and averages over the spins. 

2 ^ (v(P£+,A)|7^'|m(;?<,-,A'] 
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(34) 



where e(,e represents the coupling of the evanescent photon to the electron. Combining the spinor bilinears 



into traces, 
1 
2 
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The final result is 
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(36) 



again in agreement with Eqs. (2][3). As promised, under the DRED scheme renormalization program, 
evanescent Green functions are rendered finite by renormalization and contribute to scattering amplitudes 
at order (e). Also as promised, the results are completely equivalent to those of the CDR scheme. 



V. THE FOUR-DIMENSIONAL HELICITY SCHEME 



In the four-dimensional helicity scheme, one defines an enlarged vector space of dimensionality Dm = 
4 — 2e, in which loop momenta take values, as in the CDR scheme. In addition, one defines a still larger 
vector space, of dimensionality D., = 4, in which internal spin degrees of freedom take values. The precise 
rules for the FDH scheme are given in Ref. Q. They are: 

1. As in ordinary dimensional regularization, all momentum integrals are integrated over dimen- 
sional momenta. Metric tensors resulting from tensor integrals are D„, dimensional. 

2. All "observed" external states are taken to be four-dimensional, as are their momenta and polarization 
vectors. This facilitates the use of helicity states for observed particles. 



18 



3. All "unobserved" or internal states are treated as Ds dimensional, and the Dg dimensional vector 
space is taken to be larger than the dimensional vector space. Unobserved states include virtual 
states inside of loops, virtual states inside of trees as well as external states that have infrared sensitive 
overlaps with other external states. 

4. Both the Dy and dimensional vector spaces are larger than the standard four-dimensional space- 
time, so that contraction of four-dimensional objects with or Dg dimensional objects yields only 
four-dimensional components. 

To keep track of the many vector spaces and their overlapping domains, I give the result of the contrac- 
tions of the various metric tensors with one another. 







rg,v 


— Dm ) 




5^^ 5^,=D, = D,-D,n 






g^'^'v 








g^'SP 






= 0, 







(37) 



Like the HV scheme, the FDH scheme treats observed states as four-dimensional. In inclusive calcula- 
tions, however, where there are infrared overlaps among external states, the external states are taken to be 
Ds dimensional in the infrared regions. 

As in the DRED scheme, spin degrees of freedom take values in a vector space that is larger than that in 
which momenta take values. It would seem, therefore, that the same remarks regarding the Ward Identity 
and the conclusion that the Dx = Dg — D,„ dimensional components of the gauge fields and their couplings 
must be considered as distinct from the dimensional gauge fields and couplings would apply. That is 
not, however, how the FDH scheme is used. All field components in the Dg dimensional space are treated 
as gauge fields and no distinction is made between the couplings. It is common, however, to define an 
interpolating scheme, the "5s" scheme, in which = 4 — 2e 8r. The parameter 8r interpolates between 
the HV scheme {5r = 1) and the FDH scheme {5r = 0). Using this scheme gives one a handle on the impact 
of the evanescent degrees of freedom on the result, but not on the impact of a distinct evanescent coupling. 

It is claimed [5] that the essential difference between the FDH and DRED schemes is that in the former 
Dm > 4, while in the latter D,„ < 4. It must be this difference, then, that allows for the very different 
handling of the evanescent couplings and degrees of freedom. We shall see what impact this choice has in 
the calculation and discussion below. 
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A. Renormalization 

I will not give detailed results for the renormalization parameters of the FDH scheme. There is 
no point in doing so because, as I will show, the rules of the FDH scheme enumerated in the previ- 
ous section are not consistent with a successful renormalization program. The first sign that there is 
a problem with the renormalization program comes in the computation of the one-loop renormalization 
constants. In particular, the gluon vacuum polarization tensor splits into two independent components, 
n^" = IlAiQ^) {{-Q^g^" + Q") and = UeiQ'^) both of which are singular. This is a clear 
warning that what the FDH scheme calls the gluon is in fact two distinct sets of degrees of freedom. If I 
ignore Ilg and just renormalize ITa, I find the usual result that 

AP=^C,-iA^,. (38) 

Note that I also get this result if I take the spin average (trace) of the full vacuum polarization tensor. 
Because Ilg is weighted by a factor of 2e, its contribution to the spin average is not singular. Because the 
leading order term in the quantities being calculated is of order one, and the NLO term of order as, this 
result for the one-loop j8-function is all that is needed to compute the renormalized cross section at NNLO. 
Furthermore, the many NLO results that have been obtained using the FDH scheme have all renormalized 
using the above result for j8™^. 

When I try to proceed to the two-loop beta function, I find that both 11^ and ITg contribute singular terms 
to the spin-averaged vacuum polarization, while if I again ignore rig and renormalize Ha, I obtain the usual 
value for jSi, 

Pr'=Y4'^l-^CANf-lc,Nf. (39) 

This seems to be the choice made in Ref. (5) as they quote only the result for terms proportional to Q'^Q^, 
which would be part of my 11^. Since the standard lore has been that and af^ coincide, at least 

through second order corrections, this seems to be the most reasonable choice. Furthermore, it means that 
the conversion to will be |[5ll25l 





(40) 



As it turns out, it does not matter what choice one makes as even the one-loop result for j8™^, which seems 
safe if only because it is familiar, leads to the violation of unitarity. 
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B. Vacuum polarization in the FDH scheme 



Leaving aside the question of renormalization beyond one-loop, I will proceed with the calculation of 
the y-boson vacuum polarization. In performing calculations in the FDH scheme, it becomes apparent that 
the results are identical, term-by-term, to the calculation in the DRED scheme, except that the evanescent 
gluons are identified as gluons and the coupling is set to a^. Therefore I find that 
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Upon renormalizing such that 
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I find that 



FDH 



3 [nAmFDH] = L <! 1 + I ^ 1 jCf 

2 



1 + 



a. 



FDH 



+ 



a, 



FDH 



15 



131 11 



5 1 



/ 



1+ '^ ft-ln^ 



+ 



a, 



MS 



7t 



15 



133 11 



5 1 



/ 



+ 



a 



FDH 



^2 29 131 5 

"^^32+^^^^ 96" +^^^^12 



Cf 



1 + 




+ 



MS 



nil 

-Cll-CFCA^+CpNf— 1 ( l+3eln 



1 

12 



29 



45 






(45) 



C. Total Decay rate and annihilation cross section in the FDH scheme 



The results of the vacuum polarization calculation look to be disastrous as rig is singular at order af. 
However, the rules of the FDH scheme, enumerated above, specify that external states are taken to be 
four-dimensional. This means that the spin average of the vector polarizations is 



1 



spins ^ 



1 



Ml 



+ 



(46) 



which annihilates ^^^\pj)u- For the annihilation rate, the rules are a bit ambiguous, as they could be read to 
mean that the lepton spinors are four-dimensional but the vertex {y^) connecting them to the loop part of the 
amphtude is Ds dimensional. This would bring |^^^ into the calculation and lead to a singular result at 
order a^. However, Rule 4 could also be read to mean that the vertex sandwiched between four-dimensional 
states is also reduced to being four-dimensional. 
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Assuming this interpretation, I find that 
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(48) 



The results agree with one another, are correct through NLO and are finite through NNLO. Unfortunately, 
the NNLO terms are not correct! Because the discrepancy is finite, there remains the possibility that the 



conversion from a™^ to given in Eq. (40 1 is incorrect, although this would contradict previous re- 
sults ||5] |25l. If this were the case, then one would expect that the N^LO result would also be finite but 
incorrect. If, instead, the finite discrepancy at NNLO is the result of a failure of the renormalization pro- 
gram, the N^'LO result should be singular 



VI. PARTIAL RESULTS AT N^LO 



Although first computed some time ago, the vacuum polarization at four loops ifTSl [T4l remains a 
formidable calculation. It is only necessary, however, to look at a small part of the calculation: the terms 
proportional to the square of the number of fermion flavors, Nj. This is fortunate for a couple of reasons: 1) 
there are only three four-loop diagrams to be computed, see Fig. ([3]), (plus three more in the DRED scheme, 
where the gluons are replaced by evanescent gluons); and 2) the contributions from renormalization in the 
CDR and FDH schemes come only from the leading term in the QCD j3 -function (j3o and jSg). Thus, my 
result will not depend on how the higher order terms of the j8 -function are chosen in the FDH scheme. 



A. The CDR scheme 



In the CDR scheme, there are only three four-loop diagrams that need to be calculated. The first two are 
simply iterated-bubble diagrams and are essentially trivial. The third is slightly nontrivial, so I again use my 
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FIG. 3: Four-loop diagrams that contribute to the A^? term at N^LO. 



QGRAF-FORM-REDUZE suite of programs to address the problem. All of the four-loop master integrals 
can be found in Ref. ||26l . I find the result of the four-loop calculation to be 
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Renormalizing, I find 
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Using this term to compute the ot^A'j contribution to the decay rate and annihilation cross section as in 



Eqs. ( 14|17 1, 1 find the result expected from Eqs. ([2]j3]l. 
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B. The DRED scheme 



In the DRED scheme, there are three extra four-loop diagrams to compute, obtained by replacing gluon 
propagators with evanescent gluon propagators. I find 
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Upon renormalizing according to Eq. (24 1 and converting the coupling to a^^, I obtain 
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As for the CDR scheme, this leads to the expected result for the decay rate and annihilation cross section. 



C. The FDH scheme 



In the FDH scheme, however, I find that 



FDH 



a^Nc £ Q 

f 



An \ 



4e 



B\ 3 



1 1/23 1 \ 13453 3 53 1 / 
48^ + £l64"3^V^ W8^'"T8^'"2^'* 



nf (2) 



FDH 



a^Nc £ Q 

f 



An \ 



4e 



(53) 



25 



I renormalize according to 
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keeping only terms proportional to Nj. Such terms can only come from the j8™^ and j8™^ terms, 

so any uncertainty about /3f°^ has no effect here. The renormalized result is 
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The demand that external states be four-dimensional removes the lis term, but there is also a pole in 11^ 
and no finite renormalization to put the result in terms of can remove it. I must therefore conclude that 
the FDH scheme is not consistent with unitarity. 



VII. DISCUSSION 



In this paper, I have performed a high-order calculation in each of three regularization schemes: the 
conventional dimensional regularization (CDR) scheme; the dimensional reduction (DRED) scheme; and 
the four-dimensional helicity (FDH) scheme. Of these, the CDR scheme is by far the most widely used, and 
was, in fact, used to compute the original results that I use as my test basis. The FDH scheme has primarily 
been used to produce one-loop helicity amplitudes, although it has been used in a few cases in two-loop 
calculations and also as a supersymmetric regulator. The primary purpose of this paper was to put the FDH 
scheme to a stringent test and determine its reliability in a high-order calculation. The DRED scheme is 
primarily used as a supersymmetric regulator and is quite cumbersome for nonsupersymmetric calculations. 
It is, however, closely related to the FDH scheme and has been demonstrated HI |9l |23l |24| to be equivalent 
to the CDR scheme through four loops. A close comparison of the details of the calculations in the FDH 
and DRED schemes helps to identify where and when things go wrong with the former. 

In the cases of the CDR and DRED schemes, I have reproduced the known result for the hadronic decay 
width of a massive vector boson (or equivalently, the e^e^ annihilation rate to hadrons) through NNLO, 
and a few terms at N^LO. This represents computing the QCD corrections to the vacuum polarization of 
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the photon (V boson) through three loops, with partial results at four loops. In addition, I have reproduced 
the renormalization parameters of QCD (j3-function(s), mass anomalous dimension) through three-loop 
order. This establishes that I have theoretical control over all of the needed calculations through three-loop 
order. In order to obtain the partial N^LO result in the DRED scheme, I also needed the three-loop QCD 
corrections to the jS -function of the evanescent photon {V boson). 

The calculation of the V boson decay rate provides another instance of the equivalence the CDR and 
DRED schemes at the four-loop level |j23J. The ability to obtain the correct result using the DRED scheme 
required a delicate balance of the many extra couplings and their renormalization effects upon one another. 
Indeed, given the complexity needed to make the DRED scheme work, it seems that there should be little 
surprise that the FDH scheme, with its greater simplicity, should fail. 

Perhaps, it is worth considering how it is that the FDH scheme has been used successfully in so many 
calculations. Its most common use has been in the construction of one-loop scattering amplitudes via 
unitarity cuts, using four-dimensional helicity amplitudes as the primary building blocks. Thus, it is natural 
that it restricts observed (external) states to be four-dimensional. Because the FDH scheme defines that 
Dv > D,„ > 4, this restriction excludes evanescent fields from appearing as external states. This is very 



important because, as one can see from comparing Eqs. (30 1 and (55 1, terms involving external evanescent 
states are the most dangerous. Even though it does not renormalize evanescent states and couplings properly 
the FDH is able to get the nonevanescent part of the vacuum polarization tensor correct at NLO, while the 
evanescent part is ready to contribute a finite error at NLO. Because the DRED scheme defines 4 > D„„ the 
evanescent states are parts of the classical four-dimensional states. It would not seem natural to exclude 
them from appearing as external states. Instead, they are handled through the renormalization program 
so that their effects are removed from physical scattering amplitudes. In the FDH scheme, the evanescent 
states are instead additions to the four-dimensional states (as are the extra degrees of freedom that come 
from regularizing momentum integrals) and there is no barrier to excluding them as observed states. 

In an FDH scheme calculation, a tree-level term is strictly four-dimensional and is free from evanescent 
contributions. (Depending on interpretation, this may be a stronger condition than is given in the rules 
of Ref U], but it is the actual condition imposed if one defines the tree-level amplitude as being a four- 
dimensional helicity amplitude.) Because evanescent terms are absent at tree-level, they cannot generate 
ultraviolet poles at one loop. Even if one were to renormalize them properly, as in the DRED scheme, 
there would be nowhere to make the counter-term insertion! In fact, the one-loop contributions are not even 
finite, as the counting over the number of states (2e) makes the result of order e. This is clearly illustrated 



in Eq. (28 1. Neither a^, nor appear at LO. Therefore, the contributions at NLO are finite for and of 
order e (because of the counting over the number of states) for a^. In more complicated QCD calculations, 
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as will appear at LO and will therefore contribute an ultraviolet pole at one-loop, which will be removed 
by renormalization. a^, however, will still make its first appearance at NLO and that contribution will be of 
order e. Thus, one can expect that the FDH scheme, used as above, should be reliable for computing NLO 
corrections through finite order (£°). The error from improperly identifying evanescent quantities should be 
of order e. At NNLO and beyond however, the failure to properly identify and renormalize the evanescent 
parameters leads to incorrect results and the violation of unitarity. 

So, as suggested |15|, one of the FDH scheme's most important assets is that it defines > D,„ > 4. This 
feature is also the scheme's undoing, though not of necessity. Because the effects of external evanescent 
states can be removed (or indeed never seen) by imposing a four-dimensionality restriction, and because the 
effects of internal evanescent states therefore contribute at order e at one loop, it appears that one can simply 
ignore the distinction between gauge and evanescent terms. In contrast, because the DRED scheme must 
deal with external evanescent terms from the beginning, its advocates were forced to develop a successful 
renormalization program |l8l|9l. Extensive testing fSl |9l |23l |24l has shown that this program works to at 
least the fourth order and that it handles the effects of both internal and external evanescent contributions. 
As I remarked earlier, calculations in the DRED and FDH schemes are term-by-term identical, except for 
the identification of the couplings and propagating states. Thus, one could make the FDH scheme a unitary 
regularization scheme for nonsupersymmetric calculations by recognizing the distinction between gauge 
and evanescent terms and adopting the DRED scheme's renormalization program. This would, of course, 
do away with any notion of the FDH scheme being simple, but it would at least be correct. The FDH scheme 
would still be distinguished from the DRED scheme by the fact that > D,„ > 4, which facilitates hehcity 
amplitude calculations and, in chiral theories, improves its situation with regard to 75 and the Levi-Civita 
tensor ETl |28l . Furthermore, with a valid renormalization program, the requirement of four-dimensional 
observed states could be made optional. This would lead to two linked, slightly different, schemes, just 
like the HV and CDR schemes. This suggestion has already been made by Signer and Stockinger |[29l who 
in fact define their version of the DRED scheme to have precisely the > D„, > 4 hierarchy of the FDH 
scheme. 

Thus, in conclusion, the CDR and DRED schemes are correct and equivalent ways of performing QCD 
calculations through N^LO. The FDH scheme, however, has been shown to be incorrect and to violate 
unitarity beyond NLO when applied to nonsupersymmetric theories. It must therefore be viewed as a 
shortcut for performing NLO calculations and should only be used for such calculations with great caution. 
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Appendix A: Renormalization parameters for the CDR scheme 

To three-loop order, I find the coefficients of the j8 -function to be 

„Ms 2857 , 1415 o 205 1 , 79 , 11 , 
while the coefficients of the mass anomalous dimension are 

^ ^ 129 129 11413 _ _ 3 X / 139 3 A _ ^ 

^2 128 ^ 256 ^ 6912 ^ V64 8^7 ^ ^ V864 8^7 ^ 1728 

(A2) 

in agreement with known results ||301 - |33]| . 
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Appendix B: Renormalization parameters for the DRED scheme 

The coefficients of the QCD j3-function, P^^{af^) through three loops are: 

= I^Ca - ^Nf , = ^^Cl - ^CANf - ^-CfNf , 

„no 3115 . 1439 ,193 1 , 79 n 11 , 



where the notation is that 



,DR / ™DR \ ' / ^DR \ ^ / „DR \ /„DR\ ' /„DR' 



The last three indices of j8™,„ are omitted when they are all equal to 0. 
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The j3-function of evanescent QCD coupling, P^^{af^) is 
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The mass anomalous dimension in the DRED scheme is 
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The above results for j8°^, j8°'^ and y^^ all agree with the results of Refs. ||23ll24ll 

The QCD contributions to the j8 -function of the evanescent part of a non-QCD gauge coupling is a new 
result. I find 



^™io — Pv^oi — -Cf 



„nR 129 , 133 , /23 3 , \ , 10255 , /281 3 , \ 35 , 

„no /139 27 , \ , /331 81 , \ . 11 , / 195 27 , \ .5 

= - 3 fe) - (Z - + - + (1 - ^fe) c.cl 



oDR _ 3 oDR _ 25 _ 1 

Pye,02100 — g PVe.QKm — PVe,02001 " 



63 oDR _ 21 oDR _ 65 „DR _ ^9 

^ PVe,01101 — Y5 PVe, 01020 " PVe, 01002 ~ ~yi 

(B5) 



[1] G. 't Hooft and M. J. G. Veltman, Nucl. Phys. B44, 189 (1972). 



31 



[2] J. Collins, Renormalization (Cambridge University Press, Cambridge, England, 1984). 

[3] W. Siegel, Phys. Lett. B84, 193 (1979). 

[4] Z. Bern and D. A. Kosower, Nucl. Phys. B379, 451 (1992). 

[5] Z. Bern, A. De Freitas, L. J. Dixon, and H. L. Wong, Phys. Rev. D66, 085002 (2002), hep-ph/0202271. 

[6] R. van Damme and G. 't Hooft, Phys. Lett. B150, 133 (1985). 

[7] D. M. Capper, D. R. T. Jones, and P. van Nieuwenhuizen, Nucl. Phys. B167, 479 (1980). 

[8] I. Jack, D. R. T. Jones, and K. L. Roberts, Z. Phys. C62, 161 (1994), hep-ph/93 10301. 

[9] I. Jack, D. R. T. Jones, and K. L. Roberts, Z. Phys. C63, 151 (1994), hep-ph/9401349. 

[10] K. G. Chetyrkin, A. L. Kataev, and F. V. Tkachov, Phys. Lett. B85, 277 (1979). 

[11] M. Dine and J. R. Sapirstein, Phys. Rev. Lett. 43, 668 (1979). 

[12] W. Cehnaster and R. J. Gonsalves, Phys. Rev. D21, 3112 (1980). 

[13] S. G. Gorishnii, A. L. Kataev, and S. A. Larin, Phys. Lett. B212, 238 (1988). 

[14] S. G. Gorishnii, A. L. Kataev, and S. A. Larin, Phys. Lett. B259, 144 (1991). 

[15] P Nogueira, J. Comput. Phys. 105, 279 (1993). 

[16] J. A. M. Vermaseren (2000), Report No. NIKHEF-00-0032, math-ph/00 10025. 

[17] C. Studerus, Comput. Phys. Commun. 181, 1293 (2010), 0912.2546. 

[18] A. I. Davydychev, P Osland, and O. Tarasov, Phys.Rev. D58, 036007 (1998), hep-ph/9801380. 

[19] K. G. Chetyrkin, A. L. Kataev, and F V. Tkachov, Nucl. Phys. B174, 345 (1980). 

[20] D. I. Kazakov, Theor. Math. Phys. 58, 223 (1984). 

[21] S. G. Gorishnii, S. A. Larin, L. R. Surguladze, and F. V. Tkachov, Comput. Phys. Conmiun. 55, 381 (1989). 

[22] S. A. Larin, F. V. Tkachov, and J. A. M. Vermaseren (1991), Report No. NIKHEF-H-91-18. 

[23] R. V. Harlander, D. R. T. Jones, P Kant, L. Mihaila, and M. Steinhauser, JHEP 12, 024 (2006), hep-ph/06 10206. 

[24] R. Harlander, P Kant, L. Mihaila, and M. Steinhauser, JHEP 09, 053 (2006), hep-ph/0607240. 

[25] Z. Kunszt, A. Signer, and Z. Trocsanyi, Nucl.Phys. B411, 397 (1994), hep-ph/93 05 23 9. 

[26] P. Baikov and K. Chetyrkin, Nucl.Phys. B837, 186 (2010), in memoriam Sergei Grigorievich Gorishny, 1958- 

1988, 1004.1153. 

[27] W. Siegel, Phys.Lett. B94, 37 (1980). 

[28] D. Stockinger, JHEP 0503, 076 (2005), hep-ph/0503129. 

[29] A. Signer and D. Stockinger, Nucl. Phys. B808, 88 (2009), 0807.4424. 

[30] O. Tarasov, A. Vladimirov, and A. Zharkov, Phys.Lett. B93, 429 (1980). 

[31] S. Larin and J. Vermaseren, Phys.Lett. B303, 334 (1993), hep-ph/93 02208. 

[32] K. G. Chetyrkin, Phys. Lett. B404, 161 (1997), hep-ph/9703278. 

[33] J. A. M. Vermaseren, S. A. Larin, and T. van Ritbergen, Phys. Lett. B405, 327 (1997), hep-ph/9703284. 



